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Abstract 



The asymptotic representations of the functions Aii(x), Gi(x), Gi'(x), Ai^(x), Bi^(x) 
are obtained. As a by-product, the factorial identity (21') is found. The derivation of 
asymptotic representations of the integral dxAi{x)h{x,v) for v — > —oo and inte- 
grals, differing from it by the change of Ai(x) by Ai'(a;) or Aii(x), is presented. For the 
Airy function Ai{z), as an example, the Stokes' phenomenon is considered as a con- 
sequence of discontinuous behavior of steepest descent lines over the passes. When z 
crosses the Stokes ray, the steepest descent line over the higher pass abruptly changes 
the direction of its asymptotic approach to the steepest descent line over the lower 
pass to the direction of approach to the opposite end of this line. Therefore, when the 
integration contour, drawn along the steepest descent lines, goes over the higher pass, 
it begins or stops to go over the lower pass while z crosses the Stokes ray, and as a 
result the recessive series (contribution from the lower pass) discontinuously appears 
or disappears in the asymptotic representation of a function containing the dominant 
series. 

1 Introduction 

In the theory of quantum processes with particles in a constant or slowly varying field, the 
Airy function, related functions and integrals with these functions play an important role 
[1, 2]. In this paper [which is the Lebedev Phys. Inst. Preprint N 253 (1985)] we give the 
asymptotic expansions for some frequently occurring functions. 



2 Asymptotic behavior of Hi(2;), Gi(2;), Aii(^) and deriva- 



We define the Airy function and related functions as in [3, 4], but without the factor n 
Starting from 



tives Hi'(z), Gi\z) 




(1) 
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and integrating by parts, we have 

Hi(2) = J]g('=-i)(0)(-z)-'= + £„(z), e^{z) = i-z)-- dtq^-\t)e'\ (2) 
k=i -^0 

qit) = e-*^/^ q\t) = ~t\{t), 

q"{t) = (-2t + t')q{t), q^'\t) = (-2 + 6t' - t')q{t), .... (3) 

The relation (2) holds for any z, but it is useful only for Kez <^ —1 (more exactly 
for \z\ oo, \ph{—z)\ < 27r/3, see p. 432 in [4], when the remainder term 6n{z) is small 
compared with the terms of the sum). This is not the case for 2; » 1 in (2). So the decrease 
of terms in the sum (2) does not provide any guarantee of smallness of the remainder term 
[4, 5]. Its magnitude depends on proximity of saddle points (in general, critical points) to 
the integration path. Assuming z fixed and Re 2; -C —1 we shall increase n in (2). The terms 
of the sum first decrease in magnitude then increase, as the representation is asymptotic. It 
is clear that the sum of increasing terms, generally speaking, grows with n and this must be 
compensated by the remainder term; the representation with such a term is still exact. When 
n ^ 00 and the remainder term is dropped, we get the asymptotic series and the problem 
of interpretation of diverging part of the series arises [5]. Dingle [5] suggested that the 
remainder term can be restored just from the general term in each considered series. Indeed, 
the remainder term can be restored using, when necessary, the additional information about 
the function. For this reason it is convenient to regard the divergent sum of late terms as a 
symbolic representation of the remainder term [5]. With this understanding the asymptotic 
series uniquely represents its function and is the complete asymptotic expansion (for a certain 
phase range of z). In general, the complete asymptotic expansion consist of several series (two 
for solutions to second-order homogeneous differential equations) each of which determines 
its own function [5]. In this case the complete asymptotic expansion expresses the linear 
dependence of the initial function upon other functions each of which is represented by its 
own series in the considered sector of z. This linear dependence between the functions holds 
in all sectors of the complex variable z, but outside the considered sector the representation 
of each function by one series may be unsatisfactory, that is the remainder term may be not 
small. 

Asymptotic expansions with remainder terms exactly determine their functions. How- 
ever, even in the case when the remainder terms are represented only symbolically by the 
increasing terms of infinite series, it is still possible to require, that complete asymptotic 
expansions satisfy the same relations as the functions themselves [5] . 

Returning now to eq.(2), for n — > 00 we have [4] 

n=0 

Another and more simple way of obtaining (4) consists of expanding exp(— 1^/3) in power 
series and subsequent term by term integration. The integrable series is convergent every- 
where, but nonuniformly , so the integrated series is asymptotic. 

We consider now the real z — x ^ 1. The integrand in (1) first exponentially increases 
(up to the saddle point tc — \/x) then decreases. In accordance with this the asymptotic 
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expansion consists of a series of exponentially large terms (contribution from the saddle 
point) and, as we shall see below, of the series (4) (contribution from the lower limit of 
integration) . 

As mentioned above, complete asymptotic expansion should satisfy the same relations as 
the functions themselves. The relation 

Aii(x) = TT + m(x)Ai'(x) - m'(x)Ai(x), (5) 

where 

POO POO / 

Aii(a;) = J dtAi{t), Ai(a;) = J dt cost xt + — \, (6) 
is important for us here. It is easily verified by differentiation and using the equations 

Hi"(a;) - a;Hi(a;) = 1, Ai"(a;) - xAi{x) = 0. (7) 

We shall obtain the asymptotic expansions for Hi(a;), Hi'(a;) and Aii(a;) and check them 
up using (5) and the known asymptotic expansions for Ai(a;) and Ai{x) [3, 4]. 

The complete asymptotic expansion for Hi(x), x oo, can be obtained in the following 
manner. Substituting t = t + y/x, we have for the exponent in (1) 

^t_!!= 2^3/2 _ ^1/2^2 r = t-v^. (8) 

3 3 3 

Next we expand exp(— r'^/3) in power series and integrate term by term over r from —^/x to 
oo. If we simply extend the region of integration to the whole real axis (using the condition 
■\/x ^ 1), we get only the exponentially large terms and the relation (5) will not be satisfied. 
So we write 

/CO 3 /"O j 

rfte"*-V - / dte^'-"^, T> 1. (9) 
-T J-T 

In the first integral we use (8) and perform the indicated operations. In the second integral 
we expand exp(— t^/3) in power series and integrate term by term. Letting T oo, we get 

HiW~.VVVvf (10) 

n=0 n=0 

2 _ r(n+l)r(n+|) _ (6n-l)!! 

^ 3-^ ' 2-+in!7r 63-n!(2n - 1)!! " ^ ' 

The relation (10) agrees with the relation 

Hi(a;) = Bi(a;) - Gi(a;) (11) 

for the functions 

poo ^3 ^3 

Bi(a;) = / dt{ex.p{xt - — ) + sm{xt + -)}, (12) 
Jo 3 3 

Jroo ^3 
' dtsm{xt+-), (13) 
n 3 



for which the asymptotic expansions are known [4] 



oo 



Bi(x)~7rV^x-VV5]^, C = V^ (14) 

n=0 ^ 

Gi(a;)~y'^-^?^, x^oo. (15) 

n=0 

(Assuming Bi(a;) = |[Bi(a; + ie) + Bi(a; — ie)], e — >• 0, it is easy to show that the asymptotic 
expansion for Bi(x), x —>■ +oo, consists of only one series (14), cf. Chap. 3 Sections 3.4 and 
3.5 in [6] and p. 52 in [5].) 

The asymptotic expansions for Aii(a;) can be obtained by integration of the asymptotic 
expansion for Ai{t), which for t — > +oo has the form [3,4] 

Ai(^) - ^7rV^t-V4e-C f;(-l)"^, ( = (16) 



n=0 



Cn are given in (10'). Using the definition (6) for Aii(a;), let us replace the integration variable 
t by C- Then integrating by parts the leading term on the right-hand side of (16), we have 



/■oo /9 \ ^/"^ 



We unite the last (integral) term with the integral of term with ci on the right-hand side of 
(16) and integrate by parts again. Repeating the process, we get 

Ai,(^)~i7rV^^-3/V^^(-l)"^, C-U'/', v^oc. (17) 

n=0 ^ 

Here a„ are connected with c„ by the recurrence relation 

— fc — 

_ 41 _ 9241 _ 52 -203009 

oo-l, ^2-^77^, «3 - 210 . 37 ' ■ ■ ■ ■ (18) 

It can be shown now that the asymptotic series for Aii(a;), Hi(a;), Ai(a;), see (17), (10), 
(16) and the series for derivatives Hi'(x), Ai'(a:), 

m'{x) ~ 7rV2^v4eC y ^ + y ^^" +/|- , (19) 

S „ — n 



-.,3/2 



n=0 ' n=0 



Ai'(x) ^ -ivr^/^xV^e-^ ^(-l)"^, = "^c., (20) 

n=0 ^ 



satisfy the relation (5). (Each series of a complete asymptotic expansion corresponds to 
its own analytic function, so one can differentiate asymptotic series, see p. 21 in [4], and 
consequently the complete asymptotic expansion). To do this, we use the equality 
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Y,i-'^TCnd2s-n = 0, S > 1, (21) 

n=0 

equivalent to 

1 ^. ^(/c + |)r(/c + |)r(2. - k - i)r(2. -k + j) _ 

^' k\{2s-k)\ ^^^^ 



k=0 



which can be regarded as a result of substitution of asymptotic series for Ai(a;) and Bi(a;) in 
the Wronskian 

Bi'(a;)Ai(x) - Bi(a;)Ai'(a;) = tt. (22) 
Using (11) and (22), the relation (5) can be rewritten in the form 

Aii(a;) = Gi'(a;)Ai(a;) - Gi(a;)Ai'(a;). (23) 

This can easily be checked by differentiation and using the equations 

Gi"(a;) - a;Gi(a;) = -1, Ai"(a;) - a;Ai(a;) = 0. 

With the help of (23) and asymptotic series (16), (20), (15) for Ai(x), Ai'(x), Gi{x) and 
the series 

n=0 

obtained by differentiation, it is easy to recover the asymptotic expansions (17) for Aii(a;). 
The agreement of the series for Aii(x), obtained with the help of (23) and by direct inte- 
gration of series for Ai{x), is a check of correctness of series (15), (24) for Gi(a;) and Gi'{x). 
Now we note that the leading term of the asymptotic series for Gi'{x) given in handbook [3], 

7 

Gi'{x) ~ gg^~^ x^ oo, 

disagrees with (24) and its integration (integration of asymptotic series is always possible 
[4]) does not reproduce the leading term of series for Gi(x). 

It is worth-while to get the asymptotic expansion for Aii(x) with the remainder term. 
Assuming v > and using Ai"(a;) = a;Ai(a;), we have 

Ah{v)= x-^dAi'{x) = x-^Ai'{x)\'^ + x-^Ai'{x)dx= (25) 

J V J V 

poo pOO 

-v-^Ai\v) + / x-'^dAi{x) = -'t;-^Ai'(^;) - v-'^Ai{v) + 2 x-^Ai{x)dx. 

J V J V 
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Continuing the integration by parts with the help of the relation 

/ x-''Ai{x)dx = -^;-"-^Ai'(^;) - (n + l)v-''-^Ai{v) + 

J V 

POO 

(n + l)(n + 2) / x-''-^Ai{x)dx, (26) 

J V 

we easily find 

+ ^?-VT^ / x--'^"--^Ai(a;)cia;, (3n)!!! = 3- 6- 9- ...-3n = 3"n!. (27) 



(3n)!!! 

For fixed n and v ^ oo the remainder term is negligible compared with the terms of the 
sum. For v ^ 1 and n ^ oo we obtain again the asymptotic expansion in the form (23) 
with the series for Gi{x) and Gi'{x) given in (15) and (24). 
For V < we proceed similarly 

/CO rv rv 

dxAi{x) - / dxAi{x) = tt - / x'^dAi^x) = 
-oo J —oo J —oo 

/V 
x~^Ai{x)dx. 
-oo 

Finally we get 

A - / N Y^r (3^)' A-// ^ (3A; + 1)! 



(28) 



fe=0 



^^^-^ f x-'--'Ai{x)dx, V < 0. (29) 

Assuming f ^ — 1 and letting n — > oo, we obtain the asymptotic expansion in the form 
of the right-hand side of (5) with the series (4) for Hi(t') and the series for Hi'(t'), obtained 
from (4) by differentiation: 

Hi'(^)-yi^V^> x^-oo. (30) 

n=0 

We note that the leading term of the series for Hi'(a;), given in [3], 

Hi'(x') ~ — -a;~^, x — > — oo, 
disagrees with (30) and its integration does not provide the leading term of (4) for Hi(x). 
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The analog of (17) for v — > — oo is 

1/2 °° 

C=^(-«P, «^-oo. (31) 

Here a„ are the same as in (18). 

To get the asymptotic expansion for Gi{—x), x oo, we use (11) and the asymptotic 
expansion for Bi{x), see [3] or (35) below. Then we obtain 



rV2 



oo 



n=0 ^ ^ n=0 



c=^x^/^ (32) 



Differentiating it , we have 

00 

Gi'(-X) - 7rV2^V4 J-(_l)n[^^^^-2n3i^(^ + ^) _ ^^^^^^-2n-l ^^3^^ ^ _ 



4- 
n=0 



^3"n!(-x)3'*+2' ^ ) 



n=0 



Cn and dn are the same as in (10') and (20). 

3 Asymptotic expansions for Ai^(a:) and Bi^(x) for x 
— oo 

In this case [3] 

(-x)i/^Ai(x) ~ 5sin(C + ^)) - Ccos(C + ^), (34) 
(-x)V4Bi(x)~5cos(C + ^)+Csin(C + ^), C = li-xf', x ^ -oo, (35) 

where 

oo oo 

S = x/^5^(-l)'=C2,C-'', = x/^$](-l)'c2,+iC-''-^ (36) 

fe=0 fe=0 

So the combinations 

i[Ai2(x)±Bi^(x)]=^/;i,2(x) (37) 
are respectively the nonoscillatory and oscillatory parts of Ai^ (x) . Writing 

oo 

w,ix) = -[Ai^ix) + Bi^(.)] (Z^^ (38) 
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and substituting it in the equation 

w"'{x) - 4xw'{x) - 2w{x) = 0, (39) 
satisfied by Ai^{x) and Bi^(x) [3], we get 

(6n-l)(6n-3)(6n-5) (6n-l)!! 
^"-^ - (-1) 25n.3-n!- 

According to (34)- (36) eo — 1. Consequently 

n=0 ^ ' 

Similarly for the oscillatory part of Ai^(x) we write 

1 



oo 



W2{x) = -[Ai\x)-Bi\x)] 

92n . OA , 92n+l /• 2 



2^^(_,)3n+i/2--2C+^^^cos2C], C=-^i-^f'. (42) 

Substitution in (39) gives the recurrence relation 

(3n-5)(3n-3)(3n-l) 27n^ - 27n + 5 

- ^^-^ gn-2 + (-1) ^sTJT^gn-l. (43) 

The initial (^o = 1 and gi = are easily obtainable from (34)-(36). The rest are deter- 

mined from (43): 

5-41 _ 5- 7- 11 -59 _ 5- 7- 11- 12769 

5'2 — ~ 27 . 32 ' 3i ~ . 34 ' ~ 2^5 . 35 ' ■ ■ ■ ■ ^ 

4 Stokes' phenomenon and the choice of saddle points 

In this Section using as an example the function w{z), defined by the contour integral (45), 
we consider the Stokes' phenomenon, that is an abrupt appearance or disappearance of 
component series of a complete asymptotic expansion at certain phases of z. The complete 
asymptotic expansion consists of several series; the number of these component series is 
different in different sectors of the complex variable z. Each series corresponds to a contri- 
bution from a single saddle point and determines its own function. So in different sectors of 
the complex plane z the asymptotic expansion is determined by different number of saddle 
points. The problem of choosing of the saddle points, determining an asymptotic expansion, 
is connected with the topology of steepest decent lines and the disposition of ends of the 
integration path in the integral representation for the function. For this reasons we draw the 
integration path along the steepest descent lines and will watch over its deformation with 
the change of z in the complex plane. For definiteness we consider the solution of the Airy 
equation w" — zw given by the integral 

w(z)^i f dte-'^''+'^\ (45) 
Jc 
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where the contour C goes from infinity in the sector 7r/3 < ph t < 27r/3 to infinity in the 
sector — 7r/3 < ph t < (sectors 1 and 3 in fig. 1-6). 

We note that w{z) can be represented by the sum of two integrals over imaginary and 
real positive half- axes: 

w{z)^g{z) + f{z), g{z)^i dt e-'^''+'^\ f{z)^i dte-'^''+'^\ (46) 

Jioo Jo 

By substitution t = ir the integral g{z) is reduced to Yii{z). The integral f{z) = Gi{z) + 
iAi{z). Then from (11) it follows 

w{z) = Hi(^) + Gi(^) + iAi(^) = Bi(^) + iAi(^). (47) 

Introducing the modulus and phase of the parameter z = \z\e^'^ and the integration 
variable t — |i|e'^, we write the exponent in (45) in the form 

\ 1 
-i{zt+—) = |t|[|z|sin(^ + ^) + -|i|2sin3^] - cos((/? + e) + -|i|^cos3e]. (48) 

The steepest decent lines are determined by the constancy of imaginary part of (48), that is 
by the condition 

-\t\^z\cos{ip ^ Q) + -|tpcos3^] = Const. (49) 
3 

The saddle points ti 2 = |^|^''^exp[i(<^^7r)/2]. For the steepest decent lines, passing over the 
saddle points ^1,^2, the values of the constants on the right-hand side of (49) differ in sign: 

Const = 1^1^/2 sin ^. (50) 

Therefore the steepest decent lines from the passes ^1,^2 meet and coincide below the lower 
pass only in exceptional cases, when the constant (50) becomes zero, i.e. when phase of the 
parameter z is equal to 

^ = 0, ±|, ±1^, ... . (51) 

The lines of these exceptional values in the complex plane of the parameter z are called the 
Stokes rays. Where go the ends of steepest descent lines for |i| — > 00 ? Prom (49) it follows 

|^|cos((^ + 6') + ^|i|^cos3^ ^ 0, |t|^oo, (52) 

i.e. 

^^^oo = ±J(2A; + l), A; = 0,l,2, ... . 
6 

Of these directions those ones, for which the real part of (48) goes to —00 for \t\ — > 00 (i.e. 
sin^ < 0), just correspond to descent and not to ascent. But 

sin[±|(2A; + l)] = ±(-l)^ 
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is negative at odd k for the upper sign and at even k for the lower sign. Hence 



71 StT 



57r 



(53) 



6' 6 
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In fig. 1-6 in the complex t-planc the locations of saddle points ti,t2 and the behavior of 
steepest descent lines from the passes are qualitatively indicated for the specific values of 
phase of the parameter z, which changes in its own complex plane, see fig. 7. The arrows 
mark off the integration path in (45) going from the sector 7r/3<phi<27r/3 into the sector 
— 7r/3 < ph f < along the steepest descent fines. When the parameter z happens to be on 
a Stokes ray, the steepest descent lines, passing over two saddle points, meet one another at 
the lower saddle point. The integration path going over this saddle point either retains its 
direction or abruptly changes it by n/2. 

In the latter case the integration path passes over both saddle points; the lower one 
(situated at the break) appears at the integration path when the parameter z crosses the 
Stokes ray. This is a manifestation of the Stokes'phenomenon: when z crosses the Stokes 
ray the integration path begins or stops to pass over the second (the lower) saddle point. 

At the absence of break on the integration contour it contains only one (lower) saddle 
point. In this case both before and after z crosses the Stokes ray, the integration contour 
goes over only one (the lower) saddle point and Stokes phenomenon does not occur. 

In other words, the steepest descent lines, going over different saddle points, generally 
do not meet, but in one of the sectors their ends asymptotically approach each other. When 
the parameter z approaches the Stokes ray, these ends come closer and closer at a greater 
length. For z at the Stokes ray, they merge up to the lower saddle point; in this case the 
steepest descent line from the higher pass (active fine) approaches the steepest descent line 
from the lower pass (passive line) at the angle of 7r/2. When z goes beyond the Stokes ray, 
the active line breaks away from the half of the passive line and goes to infinity along the 
other half of the line. So, when z crosses the Stokes ray, the passive line changes limply, but 
the active one drastically reverses the direction of its approach to the end of the passive line 
in order to approach the other end of that line. If the integration path for z near the Stokes 
ray contains the active line, then at z crossing the Stokes ray the passive line is included in 
(see fig. 6, 1, 2) or excluded from (see fig. 2, 3, 4) the integration path, bringing in it or out 
of it the lower saddle point. This is the Stokes' phenomenon. 

If the integration path for z near the Stokes ray contains only the passive line of descent 
(see fig. 4, 5, 6), then at z crossing the Stokes ray, the active fine (with its higher pass) is 
not contained in the integration path, which, as before, goes only over the lower pass. 

So for z outside the Stokes ray, the integration path contains one or two lines of steepest 
descent with the ends going to infinity. Asymptotic expansion for the integral over each such 
line is obtained in usual manner by the saddle point method. For this reason, in order to 
get the asymptotic expansion for the considered integral, it is convenient to start from one 
of the expressions 



where Wij are defined by the same representation (45) as w{z) but with integration path Cij 
beginning in the i— and ending in the j— sector at infinity. These expressions hold for any z, 




(54) 
(55) 



W{z) = Wuiz) +W23{z), 
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but only one of them becomes suitable for the asymptotic expansion of the initial function 
in that sector of the complex ^— plane for which each of the contours can be drawn along 
the steepest descent hne over the single saddle point. 

As seen from fig. 1-6 for the sector < ph t < 27r/3 the contours C12, C23 can be drawn 
along the steepest descent lines over the passes t2, ti. So in this sector the formula (55) is 
suitable. For the sector 27r/3 < pht < 27r or — 47r/3 < pht < this expression is unsuitable 
as nether C13, nor C23 can be drawn along the steepest descent hne over a single saddle 
point. However, for z in this sector the contour C13 can be drawn along the steepest descent 
line over the saddle point ti. So here the expression (54) is suitable. 

If z is on the Stokes ray and the integration contour, drawn along the steepest descent 
lines, does not suffer a break at the lower pass, then such a contour goes over the one, 
lower saddle point only, see fig. 5, and for the asymptotic representation of integral w{z) the 
formula (54) is suitable. 

If z is on the Stokes ray and the integration contour, drawn along the steepest descent 
lines, suffers a break at the lower saddle point (and therefore passes the higher saddle point, 
see fig. 1 or 3), then it is difficult to obtain the asymptotic expansion for the integral over 
this path. Yet this integral can be regarded as a limit of the half-sum of integrals with 
parameters zi, Z2 lying on different sides of the Stokes ray and tending to z. Hence, for the 
asymptotic expansion we can use the formula 

w{z) = i[wi2 + W23 + W13]. (56) 

For z in the sector — 47r/3 < ph t < and \z\ ^ 1, using (54) and the saddle point 
method, we represent w{z) by the asymptotic series 

wiz) = w.siz) = Wj^'\z) + )(z) ~ W!;^'\z), (57) 

where Wn^^\z) is the sum of the first n terms of the asymptotic series w!h^\z), and Rn^\z) 
is the remainder term, small in comparison with wi^^\z) = 2Sn\z). The explicit terms of 
Sn\z) and Sn\z)^ connected by the relation Sn\z) = ±iSn\ze^'^'^^), are given in equations 
(68), (66). 

Similarly for z in the sector < ph i < 27r/3 and \z\ » 1 the function w{z) can be 
represented by two asymptotic series 

w{z) = Wi2iz) + W23iz) = (58) 

Wi^'Hz) + R^^^\z) + Wi^'\z) + R^:'\z) ^ W^^\z) + W^'\z), 

forming together with (57) the complete asymptotic expansion. 

(12) 

Near the Stokes ray ph 2; = the series Woo (z), representing the contribution from 
the higher pass, is dominant and the series Woo yz), representing the contribution from the 
lower pass, is recessive. On the ray ph ^ = 7r/3 their role in the representation for w{z) is 
equal and near the Stokes ray ph. z — 27r/3 the series w!x,'^\z) becomes dominant and the 
series W^'^\z) recessive. 

Near the Stokes ray ph 2; = the series W^^\z) in the representation (57) and the 
dominant series W^'^\z) in (58) represent the contribution from the higher pass. Hence 
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they differ only in their remainder terms, i.e. near the Stolces ray their n— th partial sums 
are equal: 

Wi''\z)^Wi''\z)^2Sj^\z), (59) 
but the remainder terms suffer a jump 

R^^'\z) - R^^'\z) = W23{z), phz = Q, (60) 

equal to the analytic function (-2). As this one is recessive, its asymptotic representation 
and the remainder term docs not suffer any jump at this Stokes ray, cf. with the asymptotic 
behavior of ^13(2;) at the Stokes ray ph z = — 27r/3, cq. (57). 

So on the Stokes ray ph. z — the asymptotic representation for w{z) is given by the 
n—th partial sum (59) of the dominant series and the remainder term 

R'n'\z) = R(^'\z) + w,,{z) = liR'^'Xz) + Rr\z) + w,,{z)]. (61) 

All three representations (61), which are the limits of representations (57) and (58) from 
each side of the Stokes ray and their half-sum, are equivalent. 

Near the Stokes ray ph 2; = 27r/3, due to dominance of the higher pass, the n—th partial 
sums Wn'^\z) and Wn^\ze~'^'^'^) are equal (the branching of the sums at 2; = is taken into 
account) : 

Wl^'\ze-'n = Wi''\z) = 2Si^\ze-'n = 2^^^^^), (62) 
and the remainder terms suffer a jump at this ray 

On the Stokes ray phz = 2it/3 the asymptotic expansion for w{z) is given by the n-th partial 
sum (62) of the dominant series and the remainder 

i?(i3)(^,-2..) ^ ^(23)(^) ^ ^ l[R^:'\ze-'n + RT\z) + w,,{z)]. (64) 

Note that the functions Wij{z) can be turned one into another by rotation of the integra- 
tion path over the angle of ±27r/3: 

w,-i{z) = -e^^'/^W2s(ze^'''^^), wu(z) = e-^^'/^W23(ze-^'''/'). (65) 

The function W23 up to a factor coincides with the Airy function: W23{z) = 2iAi(2;). As seen 
from fig. 5, 6, 1, 2, 3 for z in the sector — 27r/3 < ph 2; < 27r/3, the integration path C23, 
drawn along the steepest descent line, goes through one saddle point. So for |2;| ^ 1 the 
Airy function is represented by one asymptotic series, the recessed one near the Stokes ray 
ph 2; = 0, cf. (57). By the saddle point method we get for this series 

27r 

Ai{z)^S^^\z) + Rn{z), \phz\<—, 
^^H-) = ^e-^X:^^(-0-^ C = lz'^'. (66) 

fe=0 
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For |ph 2;| < 27r/3 the remainder term Rn{z) is small in comparison with at least first few 
terms of Sn\z). Outside this sector the value of the representation (66) diminishes as z ap- 
proaches the negative half-axis because Rn{z) increases. In order to get better representation 
for the complementary sector 27r/3 < ph 2; < 47r/3, we use the relation 

ki{z) = -e-^^'l''ki{ze-^^'''^) - e^^'''^ki{ze-^^''''^), (67) 

following from (55) and (65). For variable z in the complementary sector the arguments of 
the Airy functions in the right-hand side of (67) does not leave the main sector, where (66) 
holds. It is easy to verify that 



-2-Ki 



1/2 



fe=0 



and 
Hence 



-e2'^*/3Ai(^e-^"'/=^) = S^^\z) - e^'''I^Rn{ze-^''"''). (69) 

Ai(z) = >5«(z) + -^^(-^e-^'^O - e-2'^^/=^i?„(ze-2'^^/3) - e2"^/=^i?„(^e-^"^/=^), 

27r/3 < ph ;2 < 47r/3. (70) 

As seen from here, the remainder, represented by the last two terms in the right-hand side 
of (70), does not exceed the double maximum value of Rn{z) in the main sector. 
Similarly , we find 

w{z) = -2ie^'''/^Ai{ze^'''/^) = 2S^^\z) + If^^'I^Rnize^^l^), -47r/3 < ph z < 0, (71) 

w{z) = 2Sl^\z) + 25(2) (ze-2-*) + 2iR^{z) + 2e-''''^Rn{ze-'^''''^), < ph ^ < 27r/3. (72) 

Asymptotic representations for Bi(2;) follow from the given relations with account of 
equality ^i{z) = w{z) — iKi{z). 

As seen from (70) and (72), the asymptotic representation in the complementary sector 
is the sum of two analytical continuations of the asymptotic series from the main sector, 
continuations corresponding to the right and the left detours of the branch point z = 0. 
Thereby the asymptotic representation does not depend on the position of branch cut. 

5 Asymptotic behavior of dxAi{x)h{x,v) for v —oo 

The probabilities of many processes with elementary particles in a constant external field 
are reduced to the integral [1,2,7] 

/■oo 

/ dxAi{x)h{x,v). (73) 

J V 
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For the processes taking place also in the absence of a field, the parameter v is negative 
and V — > — oo for the field tending to zero. Let as consider the asymptotic behavior of the 
integral (73) for v — > — oo. It is assumed that the function h{x, v) has no singularity in the 
integration interval and can be expanded at x — v in the series 



h{x,v)=J2hk{x-v)'/'=J2fk(^) , (74) 

k=-l k=-l ^ ^ ^ 

in which the coefficients fk weakly depend on v, 

fk = hk- i-v)'^/' ~ /o. (75) 

In the following we omit for brevity the second argument of the function h{x,v). Then we 
can write 

poo poo 

/ dx k\{x)h{x) ^ h{Q)A.\i{v) + / dxxk\{x)h{x), 

J V J V 

= M£W^ =±hi.- vfl- = f; A (i^) . (76) 

fe=-l fe=-l ^ ^ 

The function h{x) has the property of h{x): it is finite at x = and behaves as /_i (^i^) ^'^^ 
near x = v. 

Introducing the function 

2 , ■.k/2oo , \k/2 

k=-l ^ ^ k=3 ^ ^ 

vanishing at a; = f not weaker than {x — t>)'^/^, and using equation xAi{x) = Ai"{x), let us 
transform the second term in (76) by integration by parts: 

/oo j2 
dxxAi{x)h{x) = f_,y/^—[2'/^Ai\t)] - foAi'iv)- 

_^A.i2^/^Ai'{t)]-f2^ + J^°^dxAi{x)H"{x), t^2-^/\. (78) 
Here we have used the equation [1] 



r°° dxAijx) ^ r°° dxAi{x + v) ^ 22/3Ai2(i) t = 2-2/3 

Jy y/x-V yfx 



V. (79) 



It can be shown that the coefficients are connected with in (74) by formulae 

k ^ k 



/2fe-l = -V/2.-l, f2k^-Tf2s--hiO). (80) 

V ^-^ V ^-^ V 

s=0 s=Q 

Using (80), we obtain from (76) and (78) 

/"°° ,,,,,, , r dx,., , ,Ai(t)Ai'(t) tAi''(t)+Ai'^(t), 
dxAi{x)h{x) = hm'K-2 j ^-Ai{x)]+f_^[-^^-^ 
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,Ai(v) Ai'iv), Ai(t)Ai(t) Ai(v) ,,,,,, , m 

Here hi{x) = H"{x) and in the first term the equation (28) was used. The function hi{x) 
near x — v and x — Q has the properties of h{x), but is {—vY times less than the last one: 



hi{x) - {-v)-^h{x). (82) 

Applying to the last integral on the right-hand side of (81) the same procedure as to the 
original one and repeating this process, we obtain finally the asymptotic representation 

Ai(t)Ai'(t) tAi2(0+Ai'2(t) 



rdxAi(a;)M2;) = E{/ife(0)[7r-2 f ^Ai(x)] + 



Ai(^) Ai» Ai(t)Ai-(t) Ai(t;) , _ .-v. .oo^ 

Jk,0[ 2 ^ +/fe,l— 7 W2 •/*^'2 5—), t — Z ' V. [^6 

Here, similarly to (74) we have used the notation 

m/2 



hk{x,v) = hk{x) = ^ /iik,m(2; - V)""/^ = ^ /jfc,^ ( — j 

m=— 1 m=—l ^ ' 



(83') 



/ifc(x) is obtained from hk_i[x) in the same manner as hi{x) from hQ{x) = h{x) = h{x,v). 
It is easy to obtain a connection of the coefficients fk,m with the coefficients of the previous 
function: 

_ (m + |)(m + |) -^% 

Jk,2m-1 — § Jfc-l,2s-l, 



= + A_.,. - /.._.(0)1. (84) 

With the help of these equations the odd coefficients fk,2m-i can be expressed via odd 
coefficients /2s- 1 of the initial function, and the even coefficients fk,2m - via even coefficients 
/2s of function h{x) and via values of functions hn{x), n < k, at x = 0. 

As the connection of fk,2m with /2s and hn{0) is linear, let us first find the dependence of 
fk,2m on hn{0) putting at first /2s = 0. So, using the second equation in (84) first for k — 1, 
then for A; = 2, we express /2,2m through h{0) and hi{0). Continuing this process, we find 

n=0 n=0 1=0 

, (n,-i + 2)! ^;^^ (n2 + 2)! ^;^^ (ni + 2)!V0) ^ 1^ /^.-i-.(O) 
+ 1 2^ i 1 -iuo — / — ^ Vm,s , [00) 

nfc_i=0 n2=0 ^ ni=0 ^ s=0 

_ (m + 2)! !fi^ (n, + 2)! (712 + 2)! (m + 2)! _ (m + 2 + 3^)! 

m! ^ n,! '"^ 712! ^ n,\ (3s)!!!m! ' ^ ^ 

ns=0 * n2=0 ^ ni=0 ^ ^ ^ 
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The final formula for Qm,s is easily proved by the use of known formula 

^ (n+p)\ _ (m+p+l)\ 

^ n\ m!(p+l) ^ ^ 



n=0 



and the induction method. By definition, (3n)!!! = 3 • 6 • 9 • . . . • (3n), 0!!! = 1. Thus 

fc=0 s=0 ^ n=0 

Then the combination of terms with even coefficients jk^2ra in formula (83) transforms into 
the following expression 

^ Ai(t;) Ai» ^ Ai(^)_ 

k=0 k=0 

-Y^h ^QuAi>)v (3^^ + 2)1 , Ai(^)v (3s + 4)! 

n=0 s=0 ^ ^ s=0 ^ ^ 

Using the relation Ai(x) = xAi"{x) and twice integrating by parts, we find 

y_ ^Ai(x) = ^ + (n + 1)^ + (n + l)(n + 2) |^ ^Ai(x). 

Repeatedly using this equation, we obtain 

" dx _ Ai'jv) ^ (3g + 2)! Ai(tO ^ (3g + 4)! (3n + 5)! T 

H^) - (35)11^3^+^;^^ (3s + 3)!!!^;3-^(3n + 3)!!!y_^^3;^ 



(90) 

Thus the expression in braces in (89) is the asymptotic representation of integral (90). So 
for the case f2s = 0, i.e. when only odd k are present in (74), the expression (83) becomes 
simpler {t = 2~'^/^v): 



dxAi(x)M^) = J]{7rM0) ^i;^3/;^ (/fe,-i + /M)}- (91) 

Prom the expression 



connecting the function hk+i{x) with the previous one, it follows 

^.+i(0) = lh^i\0) - + 4^(/.,-i + fk,i). (93) 
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Substituting /c + 1 by /c in (92) and differentiating three times the obtained hk{x), we find 
for X — 0: 

Continuing to use formula (92) in such a way, it can be shown that 

E (STpI-i^l^A-.-' {l)^ + ^'""'iti'""'" (f) J- (''S) 

where the products of half-intiger numbers are written in terms of Pochhammer's symbol 

N / N T(a + n) 

{a)n = a(a + l)(a + 2) ■ ■ ■ (a + n - 1) = ' , , . 

r(a) 

Then 

E = Ef^F - i A^-' ■ + =Hi^*3/.(«)}. (96) 

fe=0 ifc=0 ^ 

Thus the asymptotic representation of the initial integral takes the form 

k=0 



1 AV ^ >F^r^^^''nO) riAi2(t) + Ai'2(i) 3 , , , ^ , 



Ai{t)Ai'{t) hdv) ., 9-2/3,, (f _ni rqsi 

It only remains to express the coefficients fk,±i via the odd coefficients fis-i of function 
h{x). It can be done by repeated employment of the first formula in (84): 

(m+|)(m + |) ^ 1 3 

/fe,2m-i = 2^ (nfc_i + -)(nfe_i + -)••• 

"fc-i=0 

n3+2 ^ „ "2+2 "1+2 

• • • E + ^){n2 + 2) E + 2)^^^ + 2^ ^ ^^^^ 

n2=0 ni=0 s=0 

In this expression we reverse the order of summation first over rii and s, then over 712 and 
s, and so on. Finally, we obtain 

/ , 1\/ , 3 \ rn+2k 
/fe,2m-l = 2^ R2s-l{k,m)f2s-l, (100) 
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m+2 ^ 2 "3+2 ^ 2 "2+2 ^ 2 

nfc_i=nfc_i(s) n2=n2(s) ni=ni(s) 

(101) 

f . S 0, 0<s<2i, 
'''^'^^\s-2i, s>2i, ^2.-i(l,m) = l. 

The formula (98) is correct if the expansion of h{x) = h{x, v) at x — v contains only the 
odd coefficients /fc, see (74). Just such functions were considered in paper [2], where the 
representation (98) was obtained. For /2s 7^ it follows from the lower formula (84) that to 
the expression (85) for even coefficients fk,2m, the sum 

(102) 

analogous to the expression (100) for odd coefficients fk,2m-i, is added. In this sum 

m+2 

T2sik, m) = Yl K-i + + 2) • • • 

nfc_i=nfc_i(s) 
"3+2 n2+2 

Y (^2 + l)(n2 + 2) Y («i + l)K + 2), r2,(l,m) = l, (103) 

n2=n2is) ni=ni{s) 

and the nj(s) are the same as in (101). As a result in asymptotic representation of integral 
(83), the terms linear in Ai{v), Ai{v) remain: 

- dxMixMx) = Yi^ - [ + i^^V.(-)]/^>-+ 

MtWjt) <P3/2{v) .. . AiV) Ai(t;) f-2-y\ (^0A) 

We remind that (pa{v), a = 3/2, 5/2, is given in (97) and also that h{x) = h{x,v) and 
^(3fc)^Q^ = d^^h{x,v)/{dx)^'' l^^o is a function of w. 

6 Comments to the formula (104) 

We note now that all terms of asymptotic representation (104), except /i^^'^^(O), are oscillatory 
functions of v. Indeed, the linear in Ai(f), Ai'(f) terms for w <C — 1 arc pure oscillatory 
according to (34). The quadratic in Ai(t'), Ai'(f) terms for v <^ —1 contain both oscillatory 
and nonoscillatory parts, but the last arc canceled with nonoscillatory functions (j)a{v)- 
Really, as it was shown in Section 3, Ai^(t) is the sum 

Ai\t) ^ wi{t) + W2{t) (105) 
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of nonoscillatory wi and oscillatory W2 functions, the asymptotic series of which are given 
by formulas (41), (42). Differentiating (105), we obtain 

Ai{t)A]'{t)^^w[{t) + ^w',{t), (106) 

tAi'it) + Ai^it) = ^w'lit) + ^w'M (107) 
It follows from (41) and (97) that 

W[it) (p3/2{v) W'lit) 305/2 (t^) 



(108) 



27r(-'y)3/2 4^;3 ' 2nV^ 4^;3 ' 

so that combinations in square brackets in (104) are pure oscillatory functions: 

Ai(t)Ai^(t) 03/2 (t^) ^ w',{t) ^ 1 _ 5-7-89 

7r(-v)3/2 ^ 4^;3 27r(-t;)3/2 ^ 2^ • 3 • -u^ 2^ • 3^ • -yS ^ ^ ^ 

1 5-7 5 ■7- 11 -13 -107 , , 

+ ---]cos2C, (109) 



2(-y)3/2 26 . 33 . (-?;)9/2 ■ 3^ • (-i;) 15/2 

tAi2(t)+Ai'2(t) 3 ^ , , <(^) 



+ x::8<^5/2(^) 



TT^/^ 4^3^ ^ 27^^/^ 

r -. 1 5 •7- 3923 , . 

[-^-2^:3^- 21^3^7^ + ---]«-2C+ 

[- 23.3.U3/^ - 2^.3^'(iV ^---^"^^^' C^^(-^r4(-r/^. (110) 
Let us further simplify the coefficients i?2s-i(A;,m) and T2s{k,m.) in formulas for fk,2s-i 
and /^2s si^t least for the lower values of k, see (100-103). We begin with T2s{k,m). Using 
(87), we obtain for the first sum in (103) 

ni=ni(s) ni=0 ^ ^ ni=0 



(n2 + 2)!3 V 2y (s-3)!3' 

As s takes only integer positive values, the step function 6 may be omitted because (s — 3)! 
turns into infinity at integer negative s — 3. Putting n2 — m, we have from (111) and (103) 

_ / X (m + 5)! .s! 
For the second sum in (103) we obtain on account of (111) 

n3+2 n2+2 

J2 (^2 + l)(n2 + 2) {ni + l){ni + 2) = 

n2=n2(s) ni=ni(s) 
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K + 8)! ^!(^3 + 5)! s\ 

3- 6- (713 + 2)! 32.(s-3)!(n3 + 2)! 3-6-(s-6)!' ^ ^ 

Hence, from (103) and (113) it follows for A; = 3 

r m^ - —^!I^±31 ^'(^ + 5)! s\ 

' > 3 • 6 • (m + 2)! 32 • (s - 3)!(m + 2)! ^ 3 • 6 • (s - 6)! ■ ^ ' 

Similarly, 

(m + 11)! s!(m + 8)! 

r25(4,m) = 



3-6-9 -(771 + 2)! 32 -e- (s-3)!(777 + 2)! 
s!(7n + 5)! s\ 



(115) 



32-6- (s-6)!(777 + 2)! 3-6-9- (s-9)!' 
and so on. 

With the help of these expressions it is easy to find several first coefficients f^2m (^02), 

/+ =27;-=^(/o + /2 + /4), 

/lt2 = 6t'-=^(/0 + /2 + /4 + /6), 

/+0 - 47;-^[10(/o + /2 + h) + 9/6 + e/s], 
/2+ = 127;-6[20(/o + /2 + A) + 19/6 + 16/8 + lO/io], 
/3+ = 2^ • 57;-9[22 . 7(/o + /s + A) + 3^6 + 2=^ • 3/3 + 2 • 3^10 + 3^12], 
/3+ = 2^ • 3 • 5v-'[2^ ■ 3 • 7(/o + /2 + /4) + 2 • 41/6 + 2' ■ IQ/s + 2^10 + 3' ■ 5/i2 + 3 ■ 7fu]. (116) 
To calculate R2s-i{k,m), we proceed similarly. Using the relation [8] 

^ r{k + a) _ r(77 + a + l) Tja) 

^r{k + b) {a-b + i)r{n + b) {a-b + i)r{b-iy ^ ' 



k=0 



we obtain for the first sum in (101) 

Putting here 712 = in, we obtain i?2s-i(2, 777). To calculate i-th sum, the relation 

r(77, + 1 + 3/) _ 1 r(77,+i + I + 3/ + 3) (6/ + 3)!! 

n^his) r(n. + i) " ^(my ^ r(77.+, + 1) + ^^(^ 

/ 1 _ A r(.- 27 + 1 + 30 ^ (6[+3)!! 

% 2 r(s-27-i) + 23(w) ^-^^y^ 

is useful. 
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For the second sum in (101), we find with the help of (118), (119) 



r(n. + 1) r(H, + |) _ r(n. + f) nn, + f) i 

nhs) + i) ni:.,s) + I) 3 • 6r(n3 + f) 3r(n3 + 1)^8 



h\J{s + \) 1 9!! I nf 5V r(s + i) ^ 1 



'l^-2j [^r(^ + 8]> + ^^ + 2^-%^-2j ^ 3 -8^.-1) + 2^] + 



2^ ^3-6r(s-^) + 3.8r(s-|) 3-6-26J- ^'^^^ 

At 713 = m the right-hand side of (120) equals to it!2s-i(3, m). The first few coefficients 
fk,2s-i in (100) are 

3 15 

fo-l — f-l, f 0,1 — hi fl-l — -^{f-l + fl + f3), fl,l = -^{f-l + fl+ f3 + fb), 

/2,-i = yI;^[53(/-i + /i + /a) + 50/5 + 35/r], 
/2,i = ^[116(/_i + A + /3) + 113/5 + 98/v + 63/9]. (121) 

7 Asymptotic behavior of integrals related to the in- 
tegral (73) 

In applications the asymptotic behavior for f <^ —1 of integrals obtained from (73) by the 
replacement of Airy function Ai(x) by the function Aii(x) or Ai'(x) are needed. In the first 
case, integrating by parts, we obtain 

roo poo 

/ dx Aii{x)h{x) — / dxAi{x)g{x), h{x) = h{x,v), g{x) = g{x,v), (122) 

J V J V 

2 1 

g{x) = / dx h{x) — 2h-i{x — V {x -v)^^"^ + -h2{x -v)^ -\ . (123) 

Jv 3 2 

It is seen that g{x) has the same structure as h{x) in (74). Hence the above consideration is 
applied to the right-hand side of (122). 
In the second case we have 

POO 

/ dxAi'{x)[h_i{x - v)-^/^ + ip{x)] = 

J V 

poo 

^2h^iAi{t)Ai'{t)-hoAi{v)- dx Ai{x)(f'{x), t^2-^/^v, (124) 

J V 

h{x) = h_i{x - v)-^/2 + /io + /ii(a; - 't;)^/2 + /i2(x - -y) + • • • = h_i{x - v)-'^/^ + ip{x). (125) 
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The first term in the right-hand side of (124), arisen from the first term on the right-hand 
side of (125), is found by differentiation with respect to v of the expression 

POO POO 

/ dx ki{x){x - v)-^'^ ^ dtt-^'^ki{t + v)^2''^ki^{2-^'\). (126) 

Jv Jo 

The function (p'{x) according to (125) again has the structure of h{x) in (74) and the pre- 
ceding consideration is apphcable to (124). 
Note now that 



dxAi{x)^'{x) = -^Y.^V^ + O.T., (127) 
where O.T. are oscillatory terms. As 

V,_ (-l)"(2n-l)!! 
dx" ^ ' 2^{x - v)"+V2 ' 

then 

^(3»«)(o) = A(3»«)(o) + (_i)»_g=±ip__L_. (128) 



Hence, 



oo 



^(3n+i)(o) _ ^ /i(3"+i)(0) - (-l)»(6n + l)!!/i-i 

(3^)!!! (3n)!!! ^ ^ (3n)!!!23'^+i(-t;)3"+3/2 " ^ > 

On the other hand, for the nonoscillatory part of the term 2/i_iAi(i)Ai'(i) in the right-hand 
side of (124) we have according to (41) and (106): 



1 >^ (-l)"(6n+ l)!!/i_ 

- (t) _ TT ^3^^in23n+l(_^)3n+3/2 ' 

So, the all nonoscillatory part (130) is cancelled by the second sum on the right-hand side 
of (129). Consequently, the nonoscillatory part of (124) is given by the first sum in the 
right -hand side of (129), and the oscillatory part can be found according to the previous 
consideration. 

Note added to the electronic version of this paper. 

The consideration of the recessive series switching on and off, carried by us in the Lcbcdcv 
Phys. Inst. Preprint N 253 (1985), and connected with the abrupt behavior of the steepest 
decent line over the lower pass, led us later to the notion of the natural Stokes' line width 
[9-11]. 

The natural width of the Stokes hne is defined by the dimension of the saddle of the lower 
pass, i.e. by such change Aa of the parameter a for which the steepest decent line from the 
higher pass essentially changes its behavior near the lower pass, so that the difference of the 

phases 

Im(/i-/2) = Aa;-a + --- , Ao; = u;2(0) - u;i(0), u^,2{a) ^ -d h,2{a) / da, (131) 

of the contributions e^^^"'^ and e-'^^'^") of the lower and higher passes becomes perceptible 
quantity of the order of 1 and satisfies the uncertainty relation: 

Aa; • Act > 1. (132) 
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